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Abstract 

Using simulations to understand backgrounds from muon-induced neutrons is important in designing next-generation 
low-background underground experiments. Validation of relevant physics within the Geant4 simulation package has 
been completed by comparing to data from two recent experiments. Verification focused on the production and 
propagation of neutrons at energies important to underground experiments. Discrepancies were observed between 
experimental data and the simulation. Techniques were explored to correct for these discrepancies. 
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1. Introduction 

In ultra-low background experiments such as the 
one proposed by the Majorana collaboration [1] , it 
is critical to develop a well-understood background 
model within a verified simulation framework. Ac- 
complishing this involves developing a simulation 
package which can sufficiently simulate physics pro- 
cesses and detector response to generate an accurate 
estimation of background events for a given detec- 
tor. To this end, the MAJORANA collaboration has 
joined with the GERDA collaboration [2] to design a 
simulation and analysis software package (MaGc [3]) 
based upon the CERN software Geant4 [4,5] and 
ROOT [6]. 

Cosmogcnic backgrounds can be significantly 
reduced by performing experiments deep under- 
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ground. However, even at large depths cosmic-ray 
muons can generate untagged neutrons over a wide 
energy range (up to ~1 GeV) which may result in 
significant backgrounds. Simulations of fast neu- 
trons are difficult due to the uncertainties in and 
disagreement between theoretical models. The fact 
that primary spallation neutrons may themselves 
produce secondary neutrons and hadronic showers 
further complicates calculations and makes bench- 
marking these processes difficult. Understanding 
the background created by such neutrons involves 
the simulation of their creation, propagation and 
interaction with the detector, as well as the verifica- 
tion of all associated code. Previous works have an- 
alyzed and compared simulations with experimen- 
tal data, mainly using the FLUKA [7] simulation 
package (see [8,9]). A comparison between data and 
the FLUKA and Geant4 simulation packages has 
been completed as well [10]. Studies estimating the 
background from muon-induced neutrons have been 
performed for the CRESST experiment [11] as well 
as more generally for underground experiments [12] . 
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This work builds on previous studies (espe- 
cially [10]) to verify neutron creation and propaga- 
tion and quantify related systematic uncertainties 
within the Geant4 simulation package. New sim- 
ulations have been performed of two relevant ex- 
periments: one at CERN [13] that used a 190 GcV 
muon beam incident upon a variety of materials; 
and another at SLAC [14] that employed a 28.7 GeV 
electron beam incident upon an aluminum beam 
dump. The simulation of the CERN experiment 
is focused on validating neutron production from 
muons and the simulation of the SLAC experiment 
is aimed primarily at testing the propagation of 
neutrons and associated showers. 

Geant4 allows a user-customizable physics list 
that may be changed to suit the requirements of the 
simulation. The physics list used for the simulations 
was essentially the "QGSP_BIC_HP" list. The rele- 
vant hadronic physics utilized the Binary Cascade 
for nuclcons (below 10 GeV) and parameterized 
high- and low-energy models to describe pions. At 
higher hadron energies (>12 GeV), a quark-gluon 
string (QGS) model was employed to model inelas- 
tic interactions. High-energy and low-energy pa- 
rameterized models bridged the energy regions not 
covered by the Binary Cascade and QGS models. 
Both the Binary Cascade and QGS models employ 
a pre-compound model to describe the nuclear sys- 
tem while it is not in a state of equilibrium. At lower 
energies (0-20 MeV), neutron inelastic interac- 
tions were handled by a high-precision, data-driven 
model. The G4MuNuclearInteraction was used to 
model direct interaction between muons and nuclei. 
These interactions, also referred to as "muon spal- 
lation" interactions, involve an exchange of virtual 
photons between the muon and the nucleus and are 
a significant source of theoretical uncertainty [8]. 
Geant4 handles the virtual photons as effective pi-l- 
and pi- particles which interact with the nucleus 
using parameterized models. For each simulation, 
several versions of Geant4 from 7.0 onwards were 
used during development, but the final results pre- 
sented here are from versions 8.1 (SLAC) and 8.0 
(CERN). For more information on the various mod- 
els discussed here, please see the Geant4 online 
documentation [4]. 

2. CERN NA55 Experiment 

The NA55 experiment at CERN measured the 
double-differential cross section for neutrons cmit- 



Table 1 

Dimensions of cylindrical beam targets used in the CERN 
NA55 experiment. (See Reference [13].) 



Material 


Length 


[cm] 


Radius 


[cm] 


Density 
[g/cm 3 ] 


Graphite 


75 




i 




2.26 


Copper 


25 




5 




8.99 


Lead 


10 




10 




11.35 



ted at 45, 90 and 135 degrees from a 190 GeV 
muon beam incident upon three different materials: 
graphite, Cu, and Pb. A complete description of the 
experiment can be found in [13]. 

2.1. Simulation 

The CERN NA55 experiment has been simulated 
by Araiijo et al. using Geant4 and FLUKA and re- 
sults from these simulations are presented in [10]. 
However, that work did not consider a compari- 
son between the measured and simulated double- 
differential cross sections (i.e. cPa/dEd£l). A simu- 
lation has been completed to compare to this quan- 
tity as well as provide updated results using a recent 
version of the Geant4 package. This new simulation 
generated results consistent with those found in [10] . 

The simple simulation geometry consisted of a 
cylindrical beam target placed in the center of a 
2.2 m radius sphere filled with air. The materials and 
dimensions of the different beam targets are listed 
in Table 1. A pencil beam of 190 GeV muons was in- 
cident upon the center of the face of the cylindrical 
beam dump. For each neutron passing through the 
bounding sphere, the energy and angle relative to 
the beam were recorded. Particles were only prop- 
agated within this bounding sphere and killed once 
they exited the sphere. 

2.2. Results 

Results are presented in Figures 1, 2, and 3. 
The double-differential cross sections (d 2 a / 'dEdSY) 
were generated by selecting neutrons with energies 
greater than 10 MeV (the threshold of the detectors 
in the experiment) that satisfied 9q ± 2° for 8o equal 
to 45, 90, and 135 degrees as measured from the 
forward beam direction. The results were divided 
by the angular acceptance (2irA9), the total muons 
on target, and the cross-sectional composition of 
the material (measured in atoms/cm 2 ). To gener- 
ate the angular distribution (da/dfl), the number 
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of neutrons above 10 MeV was binned uniformly in 
cos 6 and then scaled appropriately. Our results for 
the angular dependence agree well with the results 
from both Geant4 and FLUKA reported in [10]. 

Total simulated neutron production in units of 
neutrons/muon/g/cm 2 can be determined by in- 
tegrating the results of the generated differential 
neutron-production cross section. Determining the 
total experimental neutron production is not possi- 
ble given that a very small subset of the solid angle 
space was covered. To estimate this quantity, the 
measured angular distributions were fit to a straight 
line which was then integrated. Results of this and 
comparison to simulated numbers are presented in 
Table 2. 

Since the NA55 experiment used thin targets in 
which the hadronic and electromagnetic showers 
induced by the muons do not fully develop, it is 
difficult to compare directly to results presented 
previously which considered muons through thick 
material (see, for example [8,11,10] and references 
therein). However, it is possible to compare corre- 
sponding ratios (e.g. ratio of Geant4 to experiment) 
between this work with thin targets and previous 
results with thick targets. The comparison of these 
ratios is important when trying to determine the 
source of the discrepancy between simulated results 
and experimental data. For example, equivalent ra- 
tios between simulated and experimental data for 
thin and thick targets would suggest that hadronic 
showers are handled adequately within the simula- 
tion, in [10], it is found that Geant4 is a factor of 2 
low with respect to experimental results of muon- 
induccd neutron production within liquid scintilla- 
tor (C10H20) (thick target) which has an effective 
atomic mass close to that of graphite. This ratio 
agrees well with the ratio calculated in this work 
(a factor of 2, see Table 2). For higher-Z materials 
(e.g. lead), Wulandari et al. [11] assumed that the 
neutron production versus muon energy follows the 
power law observed for liquid scintillator (~E°- 75 ) 
and then extrapolate from experimental data at 
different energies (110 GeV, 385 GeV) to an energy 
(270 GeV, LNGS depth) where simulations have 
been performed. With this extrapolation method, 
they found agreement with data within a factor of 2 
for both FLUKA and Gcant4 for lead [11]. This 
disagrees with the ratio found in this work (5.9, see 
Table 2) which suggests that Geant4 overproduces 
neutrons in hadronic and electromagnetic show- 
ers to compensate for initial underproduction of 
neutrons in lead. This conclusion is also noted by 



Araujo et al. 

The comparison between simulated and measured 
double-differential cross sections indicates a grow- 
ing disagreement as the atomic mass of the target 
material increases. In particular, simulated results 
from graphite (see Figure 1) show a reasonably good 
agreement with the energy spectrum for 45, 90 and 
135 degrees though the simulated results for 135 deg 
produce a harder spectrum. The agreement wors- 
ens for Cu and Pb (see Figures 2,3), which exhibit 
significantly over-hard spectra at 135 degrees. For 
lead and copper, the simulated angular distribution 
shows opposite curvature to the measurement. This 
was also noted in [10]. 



3. SLAC Experiment 

An experiment was performed at the Stanford 
Linear Accelerator Center (SLAC) [14] which in- 
volved a 28.7 GeV electron beam incident upon an 
aluminum beam dump. Primary neutrons were gen- 
erated mainly by photonuclcar interactions initiated 
by brcmsstrahlung photons from the decelerating 
electrons. The neutron time-of-flight and energy 
spectra were measured outside a steel-and-concrete 
shield in which the concrete width was varied. The 
final neutron spectrum outside the shield consisted 
of both primary neutrons as well as secondary neu- 
trons created through inelastic interactions of the 
primary neutrons with shield material. The exper- 
iment was simulated in detail using the FLUKA 
simulation package [15] generating results in agree- 
ment with experimental data. 

A detailed description of the SLAC beam dump 
experiment can be found in [14]. Included are 
schematic diagrams from [14] demonstrating the 
simple design employed by the experiment (see 
Figs. 4). A 28.7 GeV electron beam was incident 
upon a cylindrical aluminum beam dump inside a 
shield housing with an opening to allow the beam 
to enter. The shield housing consisted of two lateral 
steel shields and one top steel shield in addition 
to concrete positioned around all sides. An organic 
liquid scintillator detector was placed outside the 
shields laterally from the aluminum beam dump and 
the neutron energy and time-of-flight (TOF) spec- 
tra were measured for shield concrete shield widths 
of 9 ft, 11 ft, and 13 ft (2.74 m, 3.35 m, and 3.96 m). 
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Table 2 

Experimental and simulated neutron productions for the CERN NA55 experiment. Experimental values were estimated by 
integrating a linear fit to differential cross section measurements at 45, 90, and 135 degrees. 



Material 


Exp. Fluence [n//i/g/cm 2 ] 


Calc. Fluence [n//i/g/cm 2 ] 


Ratio 


Graphite 


6.8X1CT 5 


3.3X1CT 5 


2.1 


Copper 


1.2X1CT 4 


4.8xl0~ 5 


2.5 


Lead 


3.5xl(T 4 


5.9X1CT 5 


5.9 



Table 3 

A description of the composition of materials in the SLAC beam dump. Fractions are given in mass percentage and relative 
atom number as indicated. 



Material 


Density g/cm 3 


Element 


A (g/mol) 


Fraction 






O 


16.00 


50.0 (mass %) 






Si 


28.06 


20.0 






Ca 


40.08 


19.5 






Al 


26.98 


3.0 






C 


12.01 


3.0 


Concrete 


2.35 












Fc 


55.85 


1.4 






Na 


22.99 


1.0 






K 


39.10 


1.0 






H 


1.01 


0.6 






Mg 


24.31 


0.5 






Fe 


55.85 


.7 (rel. atom #) 


Steel 


7.5 


Ni 


58.69 


.2 






Cr 


52.00 


.1 


Aluminum 


2.70 


Al 


26.98 


100 % 



Table 4 

Effective A and Z for various underground sites as compiled by Ref. [12] . A comparison to concrete used in the SLAC experiment 
demonstrates its similarity in character to rock in underground sites around the world. 



Site (A) (Z) (Z)/(A) g/cm 3 



WIPP 


30.0 


14.64 


0.488 


2.3 


Soudan 


24.47 


12.15 


0.497 


2.8 


Kamioka 


22.0 


11.0 


0.5 


2.65 


Boulby 


23.6 


11.7 


0.496 


2.7 


Gran Sasso 


22.87 


11.41 


0.499 


2.71 


Sudbury 


24.77 


12.15 


0.491 


2.894 


SLAC Concrete 


24.13 


12.0 


0.497 


2.35 
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Fig. 1. Comparison between the NA55 experiment and simulations for the graphite beam target. The angular distribution plot 
(Figure 1(d)) includes data from Araujo et al.'s simulation of CERN NA55 [10] (points connected by lines) demonstrating 
that this work agrees with the previous simulation with the exception that this work exhibits an increased cross section 
at back-scattering angles. The error bars presented on the Geant4 data are statistical; error bars arc not included on the 
experimental energy spectra since these could not be obtained from a scanned image of the data from Ref. [10]. 



Concrete 




1 — 

33 ft 



Fig. 4. A beam view of the simulated geometry. The vertical 
lines in the concrete and steel materials indicate regions 
of different importance which are placed to facilitate the 
variance reduction technique importance sampling employed 
in this work. In this view, the detector is located to the right 
of the shielding. 



3.1. Simulation 

A comprehensive description of the modeled ge- 
ometry is given in [15]; the experimental geometry 
was modeled within Geant4 as closely as possible fol- 
lowing this description. A complete list of materials 
used is given in Table 3. A comparison of the com- 
position of the concrete to typical rock at selected 
underground labs is given in Table 4. 

To generate meaningful statistics within the lim- 
itations of computation time, it was necessary to 
employ importance sampling, a well-established 
variance reduction technique (e.g. see [16]). A set 
of classes implementing importance sampling is in- 
cluded within the Geant4 framework. Using this 
implementation required constructing regions in the 
lateral shield that increased in importance moving 
away from the beam dump towards the detector. A 
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Fig. 2. Comparison between the NA55 experiment and simulations for the copper beam dump. The simulated spectrum at 135 
degrees is significantly harder than the measured spectrum. The curvature in the angular distribution produced by Geant4 
exhibits curvature opposite to observations. 



neutron moving from a region of lower importance 
to one of higher importance is multiplied (i.e. the 
track is split into multiple copies) according to the 
ratio of importance values, and its "weight" is mul- 
tiplied by the inverse of this ratio. Conversely, a 
neutron moving from higher importance to lower 
importance is killed according to the importance 
value ratio. This has the advantage of effectively 
multiplying the neutron spectrum at each half 
foot of shield width and eliminating neutrons that 
travelled well away from the shield region. Geant4 
keeps track of the weight of each neutron to allow 
a comparison of the final simulation results with 
experimental data. 

Tracking cuts were established to kill particles 
cither below certain energies or in regions deemed 
sufficiently far from the lateral shield on the detec- 
tor side so as to have little or no contribution to 
the final results. In particular, cuts were established 
to kill low-energy electromagnetic particles such 



as photons (5 cm), electrons and positrons (1 cm), 
and generic charged particles (5 mm). Gcant4 used 
the distance cuts to calculate energy cuts for the 
different materials in the simulation. For example, 
these distance cuts corresponded to energy cuts 
of ~50 kcV (~4 MeV) for gammas (electrons and 
positrons) for concrete. Within steel, these cuts 
became ~200 keV (~12 MeV) respectively. In ad- 
dition, all neutrons with a kinetic energy below 
900 keV were removed. Geant4 docs not normally 
provide an energy cut for neutrons and so this cut 
was implemented within the UscrStcppingAction 
method called at each step. Since simulating the 
thcrmalization of neutrons can occupy significant 
amounts of computation time, this last cut reduced 
cpu usage. The 900 keV energy cut on neutrons was 
well below the experimental neutron energy thresh- 
old (6 MeV). This cut also eliminated gammas 
produced in neutron-capture reactions, which were 
of too low energy to contribute to the neutron flux 
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Fig. 3. Comparison between the NA55 experiment and simulations for the lead beam dump. The agreement for lead is worse 
than that for graphite and copper (see Figs. 1,2). The simulated spectrum at 135 degrees is significantly harder than the 
measured spectrum. The curvature in the angular distribution is not reproduced by Geant4. 



through photonuclear interactions. Lastly, all par- 
ticles entering a region 15 cm into the lateral steel 
shield opposite to the detector were killed, since it 
was deemed these particles would not contribute 
significantly to the flux at the detector. 

The simulations were performed for all three 
shield widths (9 ft, 11 ft, and 13 ft). As a neutron 
crossed every half foot in the shield, the energy, 
absolute time (from the initial electron event), posi- 
tion, momentum vector, and weight were recorded. 
From these parameters, we calculated total fiuence 
in addition to neutron energy spectra and timc- 
of-flight (TOF) distributions at each shield width. 
Results of these calculations are presented below. 



ered for neutrons above 6 MeV and the detector 
response was unfolded from the raw measurements 
to compare directly with calculated results. Com- 
parisons between the Geant4 calculated spectra, the 
FLUKA calculated spectra, and experimental data 
are shown in Figure 5. Geant4 demonstrates reason- 
ably good agreement with the shape of the spectra, 
but the agreement with the fiuence slightly worsens 
as the shield width is increased (this will be discussed 
further in Section 3.4). Geant4 and FLUKA both 
calculate harder spectra than the experimental data; 
this may arise due to a lack of understanding of de- 
tector response at high neutron energy and limited 
experimental statistics at these energies (see [15]). 



3.2. Energy Spectra 

The Monte Carlo energy spectra were generated 
for neutrons at or above 5 MeV, normalized per in- 
cident beam electron. Experimental data were gath- 



3.3. TOF Distribution 

Calculating the TOF distributions provided a use- 
ful consistency check for the simulation. It was found 
that the TOF distributions reproduced the shape 
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Fig. 5. A comparison of calculated and experimental energy spectra for three shield widths. FLUKA and experimental data 
were obtained from [15]. Statistical error bars are included on both of the Geant4 histograms but not on the spectra from 
FLUKA and experiment. An explanation of Reweighted Geant4 is given in Section 4.2. 



of the experimental data, but also demonstrated 
slightly diminishing flucncc with increasing shield 
width (see Section 3.4). This was consistent with re- 
sults for the energy spectra. 



taken with only the steel shield, a direct comparison 
between simulated and experimental neutron pro- 
duction is not possible. 



3.4. Fluence 

The total fluence of neutrons having kinetic en- 
ergy 6 MeV and greater was calculated for each of 
the three shield widths and compared to the val- 
ues from FLUKA and experiment (see Figure 6 and 
Table 5). Exponential fits were performed for all 
three data sets. These fits show that the attenuation 
length of neutrons simulated by Geant4 is shorter 
than measured in the experiment. The fits also in- 
dicate that the total neutron flux entering the con- 
crete shield (i.e. at a shield width of ft) is under- 
estimated by Geant4 by about a factor of 4-5. Neu- 
trons entering the concrete also pass through a steel 
shield of about a meter in width. Since no data were 



4. Estimating systematic errors of the 
simulations 

Given the inconsistencies found in the comparison 
of neutron results from Geant4 and experiment, it is 
important to quantify systematic errors and, if pos- 
sible, provide a method for correcting for these dis- 
crepancies. It should be emphasized that the tech- 
niques presented in this section are not suggestions 
for permanent fixes, but rather methods that may 
be used cautiously until as the Geant4 physics mod- 
els arc improved to better match observations. 
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Table 5 

Comparison of fluonce values for various shield widths. FLUKA and measurement results were obtained from Ref. [14]. The 
r value corresponds to the distance from the detector to the center of the beam dump. The square of this value is multiplied 
by the fluence to obtain the yield. Errors on the Geant4 fiuence are statistical. The yields for each set of results were then 
fit to an exponential using chi-square minimization to obtain the attenuation length. The results of these fits are discussed 
in more detail in section 4.2. Reference [17] cites the value of A as being 120 g/cm 2 for concrete. Reweighted (RW) Gcant4 
is a technique used to help correct for the discrepancy seen in the attenuation lengths and it is described in more detail in 
Section 4.2. The fit for RW Geant4 is further constrained by demanding that the fluences of Geant4 and RW Geant4 match 
(within errors) at zero concrete shield width. 



Source 


Shield Width [cm] 


Fluence [e ^ cm 


j 


r [cm] 


Neutron Yield [e _1 ] 


A [ff/cm^l 


FLUKA 


274 


3.59 x 10" 12 




549 


1.08 X 10 -6 






335 


9.30 x 10- 13 




610 


3.46 x 10- 7 


115 ±5 




396 


1.97 x lO" 13 




671 


8.87 x 10" 8 




Measured 


274 


2.91 x 10~ 12 




549 


8.77 x 10- 7 






335 


9.62 x lO" 13 




610 


3.58 x 10" 7 


124 ±4 




396 


1.88 x 10~ 13 




671 


8.46 x 10" 8 




Geant4 


274 


3.21 ±0.29 x 10" 


■13 


549 


9.67 ±0.88 x 10" 8 






335 


6.32 ±0.54 x 10- 


-14 


610 


2.35 ± 0.20 x 10~ 8 


97.8 ±4.4 




396 


1.14 ±0.11 x 10- 


-14 


671 


5.15 ± 0.50 x 10" 9 




Reweighted Geant4 


274 


1.07 ±0.38 x 10- 


-12 


549 


3.23 ± 1.16 x 10~ 7 






335 


1.95 ±0.77 x 10- 


13 


610 


7.24 ±2.89 x 10" 8 


119 ± 15 




396 


6.74 ±2.67 x 10- 


14 


671 


3.03 ± 1.20 x 10- 8 





4.1. Muon-induced neutron production 

A few conclusions may be drawn from the re- 
sults relevant to simulations attempting to esti- 
mate muon-induced backgrounds with Geant4. For 
forward-scattered neutrons, the shapes of the en- 
ergy and angular distributions are reasonably well 
reproduced, albeit with a production cross section 
too low by a factor of ~2 in low-Z material. For lead, 
this cross section is too low by a factor of 6. Until 
the source of this discrepancy in the physics mod- 
els can be corrected, it will be necessary to correct 
for the neutron deficit by applying a multiplicative 
factor to the normalization. The uncertainty of this 
correction should be taken to be about on the order 
of the size of the correction itself. Note that Mei 
and Hime [12] found it necessary to apply a sim- 
ilar muon-cnergy-dependent neutron-multiplicity 
correction that corresponds to a factor of 1.5 at 
190 GeV in order to correct for an apparent too-low 
spallation neutron yield for thick liquid scintillator 
targets in FLUKA. 

The situation for backscattered neutrons in 
Geant4 is much more difficult to assess, since for 
higher Z the simulated spectrum is significantly 
harder than the experimental spectrum. On one 



hand, a harder spectrum can be expected to in- 
duce more interactions with a simulated detector, 
in which case one might consider such a back- 
ground estimate to be more "conservative." On 
the other hand, for detectors that are sensitive to 
higher multiplicity events, such increased activity 
could artificially enhance the efficiency for vetoing 
such backgrounds, leading to an underestimation of 
background. Hence we urge extreme caution when 
simulations are sensitive to backscattered neutrons, 
for example in simulations of ground-shine back- 
grounds. 

4.2. Neutron propagation 

Simulation of neutron transport through a low-Z 
material (concrete) indicates a diminishing agree- 
ment with experiment as neutrons traverse a larger 
amount of material. Therefore, an estimation of 
systematic errors in neutron transport must take 
into account the dimensions of the material present 
within the simulation. Because the shape of the 
neutron energy spectrum is maintained, albeit with 
a slightly harder spectrum, it is reasonable to es- 
timate the systematic error using the integrated 
neutron fluxes listed in Table 5. For example, for 
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Fig. 6. A comparison of calculated (Geant4 and FLUKA) and measured total neutron fiuence (E> 6 MeV) for varying shield 
width. Exponential fits of the form Aexp - x / x (where x is the length through the shield) to each set of points are drawn and 
the values of the fit characteristic lengths (A) arc included in parentheses in units of g/cm 2 . See Table 5 for more details. 
Reweighted Geant4 is explained in more detail in Section 4.2. It is important to note that the technique employed in Reweightcd 
Geant4 only corrected the attenuation length (slope in the above figure) within the concrete and therefore does not correct 
for the discrepancy in total fiuence between Geant4 and data at zero shield width. 



a material of 10 ft width, Geant4 over-attenuates 
neutrons by roughly a factor of 4. Estimating this 
error for a particular experiment would involve 
some determination of average length traversed by 
neutrons, but this could prove difficult especially for 
those experiments involving non-trivial geometries. 

To the authors' best knowledge, no data is avail- 
able to which one may compare a simulated neutron 
transport through high-Z material (e.g. lead). Some 
analyses of the cascade models important for such a 
simulation have been performed [18], but these stud- 
ies have only analyzed incident protons. Therefore, 
to estimate neutron transport errors for high-Z ma- 
terial one must extrapolate from results at low-Z. 

To account for the over-attenuation of neutrons in 
low-Z material in Geant4, it is possible to "reweight" 
the neutrons as they traverse the concrete. One 
should note that this technique strays from the 
interaction model used within the Geant4 frame- 
work since the reweighting is based upon patching 
the demonstrated discrepancy between results of 



Geant4 and data. Despite this, the authors believe 
this technique provides a method by which one may 
more confidently estimate background arising from 
neutrons which must propagate through significant 
amounts of low-Z material until the source of the 
over-attenuation is addressed in the underlying 
physics models. 

The simple technique takes advantage of Geant4's 
track weighting. For each step a neutron takes 
through the concrete, the step length is used to 
calculate a reweighting factor, R: 

R = e~ ax 

where a is a constant with units of inverse length 
and x is the step length of the neutron, a is cal- 
culated using exponential fits to data and simula- 
tion results and is equal to the difference between 
the inverse attenuation lengths in each fit, i.e. a = 
1 /^exp — 1 / Ag4 is negative when the simulation over- 
attenuates neutrons. The weight of the track is mul- 
tiplied by R (which is greater than 1) at each step. 
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The weight of the neutron is propagated to any sec- 
ondary neutrons produced in subsequent interac- 
tions. This ensures that any simulation attempting 
to estimate neutron-induced background will appro- 
priately count the background from both primary 
and secondary neutron interactions. 

Results of the reweighting technique are presented 
within Figs. 5 and 6 and Table 5. Here, a value of 
a = —0.48 ± 0.09 m was used. The error on alpha 
arises from the errors on X exp and Xqa as quoted in 
Table 5. Reweighting introduces an additional vari- 
ance to the neutron flux since the technique depends 
upon the integrated path length of each neutron. 
This additional variance was estimated by analyz- 
ing the distribution of neutron path lengths and in- 
cluded in the error bars on the reweightcd total flux. 
To minimize the effect of this on the fit of the fluence, 
it was additionally required that the fluences from 
both Geant4 simulations match within errors at the 
initial concrete shield boundary. This assumption is 
justified by the fact that no reweighting was per- 
formed outside the concrete shield. The matching 
was performed using the fit to the (non- Reweightcd) 
Gcant4 results to estimate the neutron yield at zero 
shield width. This point was then added to the three 
calculated points in the Reweightcd Geant4 results, 
and the set of four points was fit to an exponential 
function. The fit for Geant4 (RW Geant4) resulted 
in a x 2 value of 0.24 (0.43) which for 1 (2) degree(s) 
of freedom yields a P-value of 0.62 (0.81). The fit 
value of A for RW Gcant4 is consistent with the mea- 
sured value within errors. 

The figures in Fig. 5 demonstrate that the 
reweighting technique effectively reduces the atten- 
uation of neutrons while maintaining the shape of 
the energy spectrum. However, this comes at a cost 
of increased variance of the overall calculated flux 
(see Figure 6) which could be partially compensated 
by increased statistics. It is important to emphasize 
that this technique could only correct for the atten- 
uation of neutrons in concrete and not for the over- 
all fluence disagreement seen at zero concrete shield 
width. Since no other measurements were made of 
the initial neutron production within the aluminum 
or neutron attenuation within the steel shield, it is 
not possible to ascertain the source of this discrep- 
ancy. However, given the underproduction of neu- 
trons generated from muons (see Section 2) and the 
over-attenuation of neutrons in concrete, it is likely 
that both contribute to the observed difference. 

Since the reweighting technique relies upon a com- 
parison of data and simulated results, it is limited 



in application to those materials for which neutron 
propagation data exists. The presented results from 
concrete may be applied to simulated rock given that 
the rock has a comparable composition to that of 
the simulated concrete (see Tables 3 and 4). 

5. Conclusion 

Gcant4 simulations involving neutrons have been 
performed and compared to experimental data. Sys- 
tematic errors within the hadronic physics models 
in Geant4 have been estimated. It has been found 
that Geant4 underproduces neutrons from 190 GeV 
muons by a factor of ~ 2 for low-Z material. The 
underproduction is more pronounced in higher-Z 
material. The shape of the muon-induced neutron 
angular energy spectra generated by Geant4 rea- 
sonably agrees with the measured shape for low- 
Z. This suggests that the underproduction may be 
corrected by applying a multiplicative factor. For 
higher-Z materials like lead, the neutron underpro- 
duction by Geant4 is more pronounced, and the sim- 
ulated backscattered neutron energy distribution is 
much harder than is observed. More work will be re- 
quired to correct the neutron production in high-Z 
materials. 

Geant4 transports neutrons adequately though 
low-Z material, though with a slightly larger atten- 
uation than is seen in experiment. This could have 
a significant effect on simulations designed to esti- 
mate neutron backgrounds within underground ex- 
periments. It is possible to correct for this over- 
attenuation by applying a reweighting technique. It 
should be stressed that this correction is intended to 
be only a temporary fix until the source of the prob- 
lems in the underlying physics models in Geant4 can 
be addressed. 
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